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Correlated photons produced by spontaneous parametric down-conversion are an essential tool for 
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quantum communication, especially suited for long-distance connections. To have a reasonable count rate 
after all the losses in the propagation and the filters needed to improve the coherence, it is convenient to 
increase the intensity of the laser that pumps the non-linear crystal. By doing so, however, the importance 



of the four-photon component of the down-converted field increases, thus degrading the quality of two- 
photon interferences. In this paper, we present an easy derivation of this nuisance valid for any form 
of entanglement generated by down-conversion, followed by a full study of the problem for time-bin 
entanglement. We find that the visibility of two-photon interferences decreases as V = 1 — 2p, where 
p is, in usual situations, the probability per pulse of creating a detectable photon pair. In particular, 
the decrease of V is independent of the coherence of the four-photon term. Thanks to the fact that p 
can be measured independently of V , the experimental verification of our prediction is provided for two 
**^ different configuration of filters. 

Pk 
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£ ; 1 Introduction 

The distribution of a pair of entangled photons to two distant partners is the building block of quantum 
communication protocols 0Gj- The entangled photons are produced by parametric down-conversion (PDC) 
in a non-linear crystal. As well-known, this process creates pairs of photons at the first order; but when the 
pumping intensity increases, four- and more-photon components become important in the down-converted 
field |s3 Ej. If one can post-select the number of photons, higher-photon components of the field may 
turn out to be a useful resource (an "entanglement laser", see 0). In other cases however, these higher- 
number components turn out to be quite a nuisance. In particular, one is often interested in two-photon 
phenomena, just think of the Bell-state measurement (BSM) that is needed in teleportation. The presence 
of higher-number components obviously degrades the quality of the two-photon interferences. In long- 
distance implementations, one can hardly overcome this nuisance by working with low pump intensities: 
after propagation along several kilometers of fibers, many photons are lost because of the losses in the fibers, 
and the efficiency of the detectors at telecom wavelengths is low, typically 10%. Moreover, if the two photons 
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come from different sources and have to interfere at a beam-splitter (as is the case for the BSM), filters must 
be introduced to ensure coherence. Thus, in order to have a reasonable count rate, one has to increase 
the pump power — and this unavoidably increases the number of unwanted higher-number components. In 
this paper, we address the degradation of the visibility of two-photon interferences due to the presence of 
four-photon components in the field, thus completing the partial study provided in Ref. 
In Section|2 we give an easy derivation of the topic and the results that is valid for any form of entanglement 
generated by down-conversion, under the assumption that the four-photon component is described by two 
independent pairs. In the rest of the paper, we relax that assumption: indeed, the four-photon coherence can 
vary from zero (two independent pairs) to one (state of single- mode down-conversion [7]) according to the 
experimental conditions [SIEl^S] • We prove that the loss of visibility does not depend on the coherence of the 
four-photon state, but only on a parameter p that is basically the probability of creating a detectable pair. 
For this full study, we shall focus on time-bin entanglement, a form of entanglement that is more robust than 
polarization for long-distance applications in optical fibers. Visibilities large enough to allow the violation 
of Bell's inequalities for two photon [HIE]; quantum cryptography jT2], and long-distance teleportation [T3| 
have been demonstrated in the recent years for this form of entanglement. For time-bin entanglement, the 
present study requires the multimodc formalism, introduced in Section [3J In fact, as the name suggests, a 
time-bin qubit is a coherent superposition of two orthogonal possibilities, the photon being at a given time 
t = (first time-bin) or at a later time t = r (second time-bin). Separate time-bins must be created by a 
pump field consisting of separate pulses: the finite temporal size (thence, the non-monochromaticity) of the 
pump pulses and the down-converted photons is a necessary feature of time-bin qubits. 
In Section we describe a setup that is used for measure the parameter p. In Section [5] we introduce 
the setup for measuring time-bin entanglement (a Franson interferometer with a suitable source) and derive 
our main prediction, namely the decrease of visibility due to the presence of four-photon components in the 
state. In Section we describe the experimental verification of our predictions. Section [7| is a conclusion. 
For readability, the technicalities of the formalism used in Sections 01 and are left for an Appendix. We note 
that the calculations of the two-photon coincidence rate provides the first explicit calculation of time-bin 
Bell experiments using the full formalism of quantum optics. 

2 Easy derivation for incoherent four-photon component 

The purpose of this Section is to derive the main results from a simple formalism, in order to gain intuition 
about the physics of the problem. The content of this Section does not apply only to time-bin entanglement, 
but to any form of entanglement obtained by down-conversion, be it with a cw or with a pulsed pump laser. 
The probabilities that we are going to introduce in this Section are "per detection window". In the case of 
a cw pump, this means "per time resolution of the detector"; in the case of a pulsed pump, this means "per 
pump pulse" ("per qubit", in the language of time-bin entanglement {14]^. 

The calculation is possible in simple terms if we neglect the coherence of the four-photon term, and assume 
that when four photons are produced, they form two independent pairs. The process of creation of indepen- 
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Figure 1: (I) Experimental setup to measure two-photon interferences, and meaningful parameters. Grey 
letters: spatial modes in the fibers; PDC: parametric down-conversions; F: filters. (11) Spectral widths of 
the down-converted photons (curves) and filters (grey shadows). Both filters are centered in the spectrum 
of the down-converted photons, and their width are < Ar- 



dent items obeys the Poissonian statistics: if P^c is the probability of creating a pair, we have P^c — \P-2c- 
For the setup, we refer to Fig. ^ We define A 0i & as the spectral width of the photons in mode a, resp. 
b, after down-conversion, that is, before the filters; the spectral width of the pump is denoted A p . As for 
the filters, we suppose that they are centered in the spectrum of the down-converted photons, and that 
they satisfy Aa,b < A Qi b to avoid trivialities; furthermore, we suppose Aa,b >> A p , so that twin photons 
certainly pass both filters, and As > Aa- Let's follow the two- and the four-photon component through the 
setup, until the coincidence detection in modes a± and 61. 

Two-photon component. To have a detection, both photons must pass the filters; because of the correlation 
in energy, if photon a passes through Fa (that happens with probability ~ Aa/ A a ), then certainly photon 
b will pass through Fb, because this filter is larger and the photons are correlated in energy. The photons 
are twins, therefore they interfere. Consequently, the detection rate due to two-photon components is (up 
to multiplicative factors) 

R2 = P 2c ^-i[f + cos(« + /3)] . (f) 

Four-photon component. Once four photons have been produced, four two-photon coincidence events are 
possible: two events in which we detect photons belonging to the same pair, and two events in which we 
detect photons belonging to different pairs. The first case is similar to the case of two-photons. In the second 
case, however, the fact that photon a passes its filter does not guarantee at all that photon b will do it as 
well; and of course, no interference will take place. All in all, 

Rl . Plt { 2 ^i [1 + C08(a+(3)]+2 |d^i^ (2) 

The total count rate is therefore R2 + Ra — R 5 [l + V cos(a + 0)] where R ~ P2c^- and where the visibility 
V is 

V = 1+ L a b =l-^+Q(P| c ). (3) 
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Recall that is singled out by the relation > . As expected, V decreases if P^c (proportional to the 
pump power) is increased. Note also that Aa < i : for a given pump power, the visibility increases if filters 
are in place. This is intuitive, considering the emission of two pairs: conditioned to the fact that photon in 
mode a has passed the filter Fa, a photon passing Fb is more likely to be its twin (whose frequency must 
lie within the filter) than an uncorrelated photon (whose frequency may lie everywhere in the spectrum). 
Finally, if only one filter is in place, then — 1 and we recover the discussion presented in Ref. [Hj. 
However, we are not really interested in fixing the pump power: rather, we'd like to fix the coincidence rate 
at the detection R. Obviously, this means that if we narrow the filters, we must increase the pump power 
in order to keep the coincidence rate constant. Strictly speaking, the quantity Pic^- is the probability per 
qubit of creating a photon pair such that the photon in mode b passes through the (larger) filter. However, 
~ holds in magnitude for typical down-conversion processes and filters; consequently, Pic^^ — R is 
an estimate of the probability of creating a detectable pair. 

The results of this Section are based on the assumption that the four-photon state is always described 
by two independent pairs. Note that this assumption is certainly good in the case of cw pump, because 
the time resolution of the detector is much larger than the coherence time of the down-converted photons. 
The assumption is more questionable in the case of a pulsed pump. The rest of the paper shows, focusing 
specifically on time-bin entanglement, that the degradation of visibility © is actually independent of the 
coherence of the four-photon term. 

3 General approach 

3.1 The state out of down-conversion 

The formalism to describe multimode down-conversion was introduced in Refs [151 for the two-photon 
component, and extended to the four-photon component for type-I down-conversion in |17| . We have applied 
this formalism to our case in Ref. [H]; we summarize here the main notations and results. 
The pump field is assumed to be classical, composed of two identical but delayed pulses: P(t) = ^JTp (p(t) + 
p(t + r)) , so in Fourier space 

P{u) = y%p{u)(l + e^) . (4) 

We use colinear type-I down-conversion in a non-degenerate regime uj s ^ tof, therefore, the signal and the 
idler photons can be coupled into different spatial modes a and b using a wavelength division multiplexer 
(WDM). The phase-matching function is written $(aj a ,cjt>); we don't need its explicit form in what follows. 
For convenience we define the following notations: 

<S>(x,y)p(x + y){l + e l{x+y)T ) = g(x, y) (1 + e^+^ r ) = G(x,y) . (5) 

The state produced by the down-conversion in the crystal reads 

|*) = iVlA^\vac) + -(A^lvac) + 0(I 3/2 ) (6) 
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where / is proportional to the intensity I p of the pump, and 

A} = J dui a du; b G(uj a , Lj b )a* (uj a )tf (u) b ) . (7) 

3.2 Detection: generalities 

We have just given the state |\E') created by down-conversion. This state evolves through the setup (in 
our case, a linear optics one so that the number of photons is conserved) according to |^) — > \ty), then 
two-photon coincidences are recorded. Here we introduce the general scheme for this detection. Let's write 
a,\ and b\ the spatial modes on which one looks for coincidences; since no ambiguity is possibly, we write a\ 
and b\ also the corresponding annihilation operators. We look at detector on mode a\ at time T A ± AT, 
where AT is the time resolution of the detectors; similarly for detection on mode b\. The coincidence rate 
reads 

rT A +AT rT B +AT 

R{T A ,T B ) = VAVB dt A dt B \\E£HtA)El+\t B )\i>)\\ 2 . (8) 

JTa-AT JTb-AT 

In this formula, 1)a,b are constant factors JS] that will be omitted in all that follows; the positive part of 
the electric field on mode a\ is defined as 

Ei+\t) = J ' dvj A (v)e-^ ai (v). (9) 

with f A (y) is a real function describing a filter in mode oi, the transmission of the filter being F A {v) = j A (y) 2 . 
The definition of E^\t) is analogous. We choose the origin of times in order to remove the free propagation 
from the crystal to the detectors. Therefore, the first time-bin at the detection is given by tj = 0, the second 
time-bin by tj — r and so on. 

Actually, formula JSJ for detection is exact for proportional counters, in which the probability of detection is 
the intensity of the field. For photon counting with a detector of quantum efficiency ry, the probability of the 
detector firing, given that n photons imping on it, is not nr\ (proportional to the intensity) but (1 — (1 — rj) n ). 
Now, for the wavelengths that we consider, the quantum efficiency is r\ ~ 0.1; moreover, the mean number 
of photons that imping on a detector is much smaller than 1 because of the losses in the fibers and in the 
coupling; finally, in our formalism we restrict to the four-photon term, so that at most two photons can 
imping on the detector. All in all, the approximation (1 — (1 — ij) 2 ) ~ 2ry holds and we can indeed use JHJ 
to compute the coincidence rate. 

3.3 Important parameters 

As we said in the introduction, we shall postpone the detailed calculations to the Appendix. All the results 
of the Sections 01 and El can be formulated using the following parameters: writing dw = duj a du>b, 

J = J duj_\g(iv a ,uJb)\ 2 , (10) 

J A = J duj_F A (u; a ) \g(uj a ,uj b )\ 2 (11) 

Jb = J duj_F B (uj b ) \g(uj a ,uj b )\ 2 (12) 
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Figure 2: Schematic of the setup used to measure the parameter p. 




duj_F A (u} a )F B (uj b ) \g(u a ,u>b)\ 



(13) 



J 4 = 



dudu/ F A (uj a )F B (Lu b ) [g*(uj a ,uj b )g*{u' a ,Lo' b )g(u a ,w' b )g(uj' a ,u} b ) + c.c] , 



(14) 



The first four numbers can be given an intuitive meaning. In fact, up to multiplicative factors: J is the 
probability of producing two photons in one pump pulse, irrespective of whether they will pass the filter or 
not; J a and Jb are the probabilities of producing two photons in one pump pulse, and that the photon in 
mode a (resp. b) passes through the filter; Jab is the probability of producing two photons in one pump 
pulse and both photons pass the filter. The interpretation of J4 is somehow more involved: it is a coherence 
term, due to the fact that the four photon state cannot be described as two independent pairs 9 . 
Obviously, Jab — J A if no filter is applied on B. But Jab = J A holds to a very good approximation also if 
A a < As, where Ax is the width of filter Fx, provided that both filters are larger than the spectral width 
of the pump A p (as we supposed in Section |2 and as will be the case in the experiment). In fact, in this 
case, detection of a photon in filter A automatically ensures that its twin photon has a frequency within the 
range of filter B, which means Fb(uji,) = 1 for all u b compatible with the phase-matching condition. 

4 A calibration setup 

Before describing the measurement of two-photon interferences (next Section) we present an experimental 
setup that allows to measure the probability p of creating a detectable pair in a simple way. This setup 
(see Fig. |2} has been presented in detail in section IV of Ref. [Sj- We give here a brief description. A 
Fourier-transform-limited pulsed laser is used to create non-degenerate photon pairs at telecommunication 
wavelengths (1310 and 1550 nm) by parametric down-conversion in a non linear-crystal. The two photons 
are separated deterministically using a wavelength-division multiplexer (WDM) and each photon is detected 
by single-photon counters (avalanche photodiodes). The signal from the two detectors are then sent to a 
Time-to-Digital converter, which is used to determine the histogram of the differences in the time of arrival 
of the twin photons. 
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We apply our formalism to this setup. For the detection, since there is no evolution but the free propagation, 
we have simply a\ — a and b\ = b. For the preparation, at first sight it seems that our formalism should 
be modified: we are dealing with a train of N pulses instead of only two pulses, so (f + e MT ) should be 
replaced with "Y^^q e lukT in formula |@}. However, a closer look shows that we can do the calculation 
without any change. In fact, in this particular setup there is no interference: then, Rc is simply the sum 
of the coincidence rates obtained when the two photons arrive at the same time, while Rl is the sum of 
the coincidence rates obtained when photon in mode a arrives a time r later than the photon in mode b. 
Since moreover R(kr, fcr) = R(0,0) and R((k + l)r, fcr) = R(r,0) for all k, we obtain R c = N R(0,Q) and 
Rc = (N — 1) R(t, 0), so we can focus on only two successive pulses. By the way, i?(0, 0) is proportional to 
the probability per pulse of creating one detectable pair (a pair that will pass the filters). 
The calculation is given in the Appendix, and the results are i?(0, 0) = I Jab+0(I 2 ) and R(t, 0) = I 2 JaJb, 
that are indeed what one expects because of the meaning of the J's f subsection 13. Therefore, in the limit 
of large N, the ratio p between the integrals of the side peak and the central peak is 

P = J -f • ( 15 ) 

Jab 

In most cases, p has a simple interpretation. In fact, whenever condition A p << < holds, we have 
seen above that Jab = J A and consequently p = IJb is the probability per pulse of creating a pair such 
that the photon that meets the largest filter will pass it. In particular, if there is no filter on mode b, p is 
the probability per pulse of creating a pair, as noticed in the Appendix of That derivation shares with 
the present one the hypothesis of small detector efficiency, but is otherwise rather different: in our previous 
paper, we supposed that a 2A r -photon state is actually N independent pairs; here, we limit ourselves to 2 
and 4 photons, but derive the result without any assumption about the coherence of the 4-photon state. 
Moreover, as argued in Sectional since J a — Jb normally holds, at least in magnitude, then p ~ Ua = 
R(0, 0) is an estimate of the probability per pulse of creating a detectable pair. 



5 The Franson interferometer 
5.1 Description of the setup 

We turn now to the main setup, which is the interferometer that allows the analysis of time-bin entanglement 
(see FigOU. This is essentially the interferometer proposed by Franson to study energy-time entanglement 
19 , completed with an unbalanced interferometer before the crystal (the pump interferometer). A laser 
pulse is first split in two in this interferometer. At its exit, we have two laser pulses with a fixed phase 
difference separated by a time t corresponding to the path length difference between the long and the 
short arm of the interferometer. In the non linear crystal, we therefore create a photon pair in a coherent 
superposition of two time-bins. 

After the crystal, the photons are separated with the WDM and each sent to a fiber interferometer in order 
to make a two photon interference experiment. 
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Figure 3: Schematic of the setup used to measure two photon quantum interference with time-bin entangled 
qubits. In addition to the two-photon coincidence, a coincidence with the pump laser provides the origin of 
time needed to define the three time-bins. 



5.2 Evolution 

The evolution of modes a and b in each arm of the interferometer is given by the following expressions: 

a f H — ► at( w ) = S{u,a)a\(uj) + C(uj,a)al(Lu) (16) 
b\co) — > gt( w ) = S{u>,p)b\{u) + C{u>,P)b\{u>) (17) 

with H3 

i _ i(ujT+ff) i I i(ujT+6) 

S(u; 7 0) = , C(u,,6) = l — — . (18) 

The evolved state \iff) is obtained by inserting the evolved operators a' and w into \^>). 
In these formulae, we have already supposed that the analyzing interferometers are identical to the pump 
interferometer. Thus, three time-bins are defined by the setup. The first time-bin, t = 0, corresponds to the 
time of arrival of photons produced by the first pump pulse and not delayed. The second or intermediate time- 
bin, t = r, corresponds to the time of arrival, either of photons produced by the first pump pulse and delayed, 
or of photons produced by the second pump pulse and not delayed. The third time-bin, t = 2r, corresponds 
to the time of arrival of photons produced by the second pump pulse and delayed. Interferences will only 
be seen when both photons arrive at t = r, because only in this case two indistinguishable alternatives are 
available. 

5.3 Two-photon interferences 

We study the detection for modes ai and b\\ all the other cases can be treated in the same way. The 
coincidence rate R{Ta,Tb) is the sum of two terms corresponding respectively to the two-photon and the 
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four-photon terms: 

R2(T A ,T B ) = I J dt A dt B \\Ei + i \t A )E { +\t B )A ] \vac)\\\ (19) 

R±{T A ,T B ) = T - J dt A dt B \\Ei+\t A )Ei+\t B )(A f ) 2 \vac)\\ 2 . (20) 

Indeed, the two- and the four-photon states do not interfere (in principle, one could insert a non-destructive 
measurement of the number of photons just after the crystal, and this would not modify the rest of the 
experiment). 

The calculation is presented in the Appendix. As said above, interferences will appear only in the interme- 
diate time-bin T A = T B = r, in which case one finds |21j : 

R 2 (t,t) = I J AB (l + cos(a + /?)) , (21) 
R 4 (t, t) = I 2 \(2J AB J + J 4 ) (1 + cos(a + /?)) + 2J A J B ] . (22) 

The result for R.2{t,t) is the expected one: one pair is produced, it passes the filters, and since it is in a 
superposition of being in both pulses it gives rise to full-visibility interferences. In the formula for R±(t,t), 
two contributions are also expected from the intuitive view of the four-photon state as two independent pairs: 
(i) the term containing J AB J means that two pairs are created, the photons of the same pair are detected 
and therefore one has full visibility; (ii) the term containing J A J B means that two pairs are created, the 
photons of the different pair are detected and therefore they don't show any interference. The remarkable 
feature is the position of the correction due to the coherence in the four-photon term, J4: it contributes to 
a full-visibility interference as well. This couldn't have been guessed without the full calculation. 
Summing 1)21(1 and l)22[l. the total two-photon coincidence rate in the intermediate time-bin reads 

R{t,t) = R 2 {t,t) + R a (t,t) = R[l + Vcos(a + [3)] (23) 

where the average count rate R is given by |U 

R = U AB + 0(I 2 ) (24) 

and the visibility V is given by V — i+i^j+^j^j^igf/jj^g) ■ Now, the terms 0{I 2 ) in the visibility are 
meaningless, because the six-photon term that we neglected completely contributes to the same order; so we 
have to keep only the first-order development of V in /, that leads to the remarkable relation 

V ~ 1 _ / 2 J±ll± = 1 _ 2p (25) 
Jab 

where p is exactly the same as defined in (I15f) . As announced in Section J4 drops out of the visibility: to 
the leading order in J, the loss of visibility is independent of the coherence of the four-photon state. 
Since p is basically the probability per pulse (so that 2p is the probability per qubit |14j ) of creating a 
detectable pair, it defines the detection rate up to multiplicative factors. Relation Q25JI therefore says that, if 
we fix a detection rate, we shall find a given visibility, no matter whether the rate was obtained by pumping 
weakly and putting no filters, or by pumping strongly and putting narrow filters. This is a positive feature: 
filters, while being useful to improve the coherence whenever this is required, do not degrade the visibility. 
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As described in sectional p can be measured independently, the relation (|25|l can be experimentally tested. 
This is the object of the next Section. 

6 Experimental verification 

In this section, we present an experimental verification of Eq. Q25JI. Two- photon interference fringes are 
recorded for different value of p, corresponding to different values of pump power, with the Franson setup 
described in the previous Section. Let us remind the reader that the down-converted photons are at the 
two telecom wavelengths, 1310 nm and 1550 nm. The measurement are reported for two different filters 
configurations; in both cases, the larger filter is on the photons at 1550nm, so this is "mode 6". 
In the first configuration, only the photon at 1310 nm is filtered with 40 nm FWHM. These data are taken 
from 0. In the second configuration, both photon are filtered. The photon at 1310 nm is filtered with 10 nm 
FWHM, while the photon at 1550 nm is filtered with 18 nm FWHM. The coefficient p is measured using the 
side peaks method explained in Section 0] The visibility for the two experimental configurations is plotted 
as a function of 2p in Fig. ^ The error bars on the experimental points represent the accuracy of the fit of 
the recorded interference patterns with a sine law J23 • The two solid lines are straight lines with slope — 1, 
according to Eq. (|25|l : the small shift between the two curves is due to the fact that the maximal visibility 
was not the same for both experiments and was left free as a fitting parameter. We observe a good agreement 
between theory and experiment. These results confirm that the loss of visibility due to four- photon events is 
directly related to p, regardless of the filtering that is applied on the photons and regardless of the coherence 
of the four-photon component This is therefore a general result very useful to estimate the effect of 
multi-pair creation in an experiment in a very simple way. 

7 Conclusion 

In summary, we have found a quantitative prediction for the loss of two-photon interference visibility due to 
the presence of a four-photon component in the down-converted field. The loss of visibility l|25(l is determined 
by the parameter p (JTSJ, that is close to the probability of creating a detectable pair. This parameter can 
be measured independently, thus allowing a direct experimental verification of our prediction. While the full 
calculation was worked out for time- bin entanglement, we have presented in Section 0a simplified derivation 
that gives the same result and applies to any form of entanglement generated by down-conversion. 
We acknowledge fruitful discussions with Antonio Acm, Christoph Simon and Wolfgang Tittel. 
Note added in proof. Since this work was finished, we have learnt of two independent papers |23l I24| that 
discussed the loss of visibility of two-photon interferences due to the presence of higher-photon-number 
components. Both calculations concern entanglement in polarization and have been done in the single-mode 
formalism: this allows to take into account the contribution of all more-photon terms and not only of the 
four-photon one. The results are compatible with ours in the regime where they can be compared (small 
pump power). Consider for instance Ref. |24j : from their eq. (2), we see that the probability per qubit of 
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Figure 4: Visibility as a function of 2p for two different filtering configurations. Full squares are experimental 
points with a 40 nm filter at 1310 nm and no filter at 1550 nm (data taken from 6 ). Open circles are 
experimental points with a 10 nm filter at 1310 nm and a 18 nm filter at 1550 nm. The solid curves are 
straight line with a slope —1, according to Eq. I|25|> . 



producing a pair is 2tanh 2 rj cosh 4 r w 2r 2 . Then our formula © predicts V ~ 1 — 2r 2 for small values of 
t, which indeed fits correctly the curve of Fig. 5 of Ref. [22] up to r « 0.5. 

A Appendix 

We recall the definitions of i?2 and R4, formulae i|19|) and l|20|) . Although we introduced them only in Section 
|5j the same quantities can be defined for the setup of Section and in fact for any setup: a given setup 
will be characterized by the relation between the preparation modes a, 6, and the detection modes ax, bx, a 
relation encoded in the operator . In this appendix, we start by working out more explicitly the general 
formulae for Ri and R4; we subsequently describe the strategy that allows a simplification of these formulae 
(a strategy already introduced in Ref. (Sj), and finally compute the explicit results announced in Sections 0] 
and El 

A.l General formula for R 2 

The calculation of R2 goes as follows. The commutation rules between the input modes a, b and the detected 
modes 01, b\ read 

ax{v a )a f {u) a ) = a) (u) a )ax(v a ) + E(uj a , a) 6{u> a - u a )l (26) 
bx(v b )P(u> b ) = P(u b )bx(v b ) + E(uj b) {3)8(uj b -v b )l (27) 

where E{uj,^) is a function that depends on the evolution undergone by the modes from the production to 
the detection — specifically, E(uj,j) = 1 for the calibration setup, while E(uj,-f) = S(uj,j) given in fll^l for 
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the Franson setup. Using these commutation relations, one finds immediately E { ^\t A ) E { +\t B ) $\ vac) = 
c(t A , tB)\vac) where we have introduced the complex number 

c{t A ,t B ) = I d^f A (Lu a )f B (Lo b )G(Lu a ,Lu b )E(Lo a ,a)E(Lo b ,l3)e-^ tA+ ^ ts K (28) 



Consequently, R2 — I J dt A dt B \c{t Al t B )\ 2 ■ The integration over t A and ts can be performed before the 
integrals over the frequencies, so finally 

R 2 (Ta,Tb) = ljdvdu/ f A (u a )f B (u) b )g{u a ,ui b )f A {u)' a )f B (J b )g*(J a ,u' b ) 

Xe -i( Ua -u' a )T A e -i{u> b -u>' b )T B ( AT )2 sinc [( Wa „ W ;)AT] sinc[(w fc - w£)AT](29) 
where we have defined the shortcut 

£(w a ,u b ,u' a ,uj b ) = E(u a ,a)E(bJ b ,p)E*(u' a ,a)E*(uj b ,f3). (30) 

In 129J1 . the first two lines are simply the expansion of the G's and the evolution term £; the last line is the 
result of the integration over t A and t B . 

A. 2 General formula for i? 4 

The calculation of R4 follows exactly the same structure as the calculation of R2, only the formulae are 
heavier. Using 

ai(v)tf(uj)tf(cj') = a^u)a^(uj')ai(u)+E(u),a)S(oj-v)^(u)') + E(uj',a)8(uj' -v)tf{uj) (31) 
and the analogous relation for mode b, one finds 

E£Kt A )El+\t B )(Xf\vac) = \AB) + \A'B') + \A'B) + \AB')=2(\AB) + \AB')). (32) 
We have defined 

\AB) = J dudu/ G{uj a ,uj b )G(uj' a ,uj' b ) Z(u a ,u b ,uj' a ,cJ b )\vac) 

where Z is the non-normalized two-photon creation operator 

Z(w„, w 6> u' a ,w' b ) = /aK)/bK) e- l ^ tA+ ^E{u: a , a)E(u> b , 0) a{ {w' a )b\(u' b ) ; 

\AB') is obtained from \AB) by replacing uj b <-> Lu' b in Z, or equivalently by relabelling the integration 
variables: 

\AB') = J (kudu/ G(uj a ,uj b )G(cj' a ,uj b ) Z((j a ,u} b ,u' a ,u' b )\vac) . 

Obviously, by simply exchanging primed and unprimed integration variables, \AB) = \A' B') and \AB') = 
\A'B), whence the r.h.s. of (|32(l . Inserting l|32|l into (|2C)|> . we see that the quantity that we must compute is 

R 4 {T A ,T B ) = I 2 J dt A dt B [(AB\AB) + {{AB\AB') + c.c.) + (AB'\AB')] 

= R4 A (T A ,T B ) + [RaATa,T b ) + c.c] + R i<3 (T A ,T B ) . (33) 
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The integrals over t A and t B are exactly the same ones that we had in the calculation of R 2 . 
The first term of the sum (|33[) is the easiest one, and can be given in closed form. In fact, in \AB) the integrals 
over uj_ and u/ are factored, and in addition, the integral on lo gives the same c(t A ,t B ) that we met in the 
calculation of R 2 , formula (|2"%|l . That is, \AB) = c(t A ,t B ) J du/G(ui' a , uj' b )a\{uj' a )b\{uj b )\vac) . Consequently, 
R 4A (T A ,T B ) = IR 2 {T A ,T B ) J du\G{oj a ,Lu b )\ 2 where we recall that G(w a ,uj b ) = g(w a , w 6 )(l + e i ^+^ T ). 
Anticipating over the discussion of the next subsection, we use here the fact that the terms that fluctuate 
in r average to zero; so the last integral is finally equal to 2 J duj_ \g(u> a , o-i b )\ 2 — 2 J. In conclusion, the first 
term of the sum (133(1 is 

Ra,i(Ta,T b ) = 2IJR 2 (T A ,T B ). (34) 
The second term of the sum J33J), R 4i2 (T A ,T B ) = I 2 J dt A dt B (AB\AB') gives 

I 2 J dudu/dQ f A (u a )f A (uJa)fB(^b)fB{^b)9*(^a,^b)9*{^' a ,^b)9(^a, uj b )g(ui a ,u} b ) 

x _|_ e -i(iia+w 6 )r~j ^ _|_ e ~i(uj' a +uj' b )T^j ^ _|_ e i(w a +W6)rj / j -f e i(«>a+&b)T^ 

x£(w a ,u} b ,uj' a ,uj b ) 

Xe -<(«.-«4> e -i( Wb -wi)T ( AT )2 sinc [( WQ _ c4)AT] sinc[(w b - wJ)AT] . (35) 
The third term R 4 , 3 (T A ,T B ) — I 2 J dt A dt B (AB'\AB') gives 

I 2 J dwdu/dti f A {u a )f A {u' a )f B {u b )f B (J b )g%Cu a ,J b )g*(w' a ^ 

x £(cj Q ,w fc ,^,w b ) 

xe _i( Wa -^)r e -i( W6 -«J)T ( AT )2 sinc [( Wa _ ^)AT] sinc[(w b - u;' b )AT} . (36) 

Note that, as it should, the difference between R^ 2 (T A ,T B ) and Ra^{T a ,T b ) is only in the contribution of 
the G's, lines one and two. 

A. 3 Strategy of the calculation 

One cannot go beyond the formulae that we just derived for R 2 and R4 without specifying what the evolution 
S(uj a , uj b , io' a , u)' b ) is (that is, without specifying the setup) and without a simplification strategy. Here is how 
this strategy goes [§]. 

1. We first notice that the times T A , T B of interest are typically 0, r etc; and as we said, £(ui a , uj b , uj' a , uj' b ) 
is also a product of terms containing either 1 or some e MT . So all our integrals for R 2 and R4 are in fact 
sums of integrals of the form J duj_J-(uA e jr2 (— ' T . Here, J 7 is a product of g(uj, o/)'s (spectral function 
of the pump and phase matching conditions) and of cardinal sines sinc((cj — lu')AT) associated to the 
time-resolution of the detectors; Q is an algebraic sum of the some of the integration variables lo_. 

2. Because we want the two pump pulses to be well-separated (well-defined time-bins), it turns out that 
all the integrals in which 1] ^ will average to zero. In fact, the typical width of g is tP } mp > 
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where t pump and At are, respectively, the coherence time and the temporal width of each pump pulse 
p(t). If the time-bins are to be well-defined, we must impose r >> At. Moreover, if one wants to 
distinguish the time-bins at detection, one must also have a sufficiently small time-resolution for the 
detector; so t » AT. In summary: in the frequency domain (which is the integration domain), if 
Q =/= the term e in ^ T fluctuates with period while in this range J-(ui) is almost constant. The 
second step of the calculation consists then in going through the factors to sort out those integrals in 
which £1 = 0. This is the clever trick that allows one to obtain readable formulae. 

3. This being done, one can also perform the limit AT — ► oo, leading to sinc(a;AT) ~ -Ap6(x). In fact, x 
is of the form ui — a/, and this is in average close to the spectral width of each down-converted photon 
But a detector cannot detect a photon unless AT >> t v c h . This is the precise meaning of the 
formal limit AT — ► oo. Obviously this limit must be performed after the estimate described in point 
2. 

In summary, for each of the setups that we want to study, we must replace £(u> a , LOb,uj' a ,uj' b ) with its explicit 
value, then by inspection identify those terms for which the dependance in r identically vanishes under the 
integral. 

A. 4 Calculations for section |3| 

For the calibration setup of section^ we must compute R(0, 0) and R(t, 0). Here, E (u> a , tot,, ui' a , u>' b ) = 1 
because the modes don't evolve from the preparation to the detection. Then, i?2 given in J2H is the sum of 
four integrals because of the product in line two; while R^2 and given respectively by (|35(l and (|36|l are 
the sum of sixteen integrals because of the products in line two. 

Let's set Ta = Tg = 0, and look first at i?2- Only the product lxl gives an integral whose argument 
does not contain r, so we forget about the three other integrals. Through the limit sinc(xAT) ~ -^S(x), 
we obtain ljj — ui'j and consequently i?2(0, 0) = I Jab- So we have obtained R(0, 0) = I Jab + 0(I 2 ) as 
announced. 

If we set Ta = r, Tb = 0, it is easy to become convinced that none of the four integrals that compose R2 
can become independent of r; therefore, R2{t, 0) = and we must compute R{t, 0) = Ra{t, 0). Obviously, 
^4,1 (t")0) = because it is proportional to i?2- By inspection, one sees that -££4,2(1", 0) = as well since 
none of the sixteen integrals can be made independent of r. In i?4 j 3(r, 0), only the integral associated to the 
product lele (with obvious notations) is independent of r. For this integral, we can again set ojj — uj'j and 
we obtain finally R(t, 0) = i?4,3(r, 0) = I 2 JaJb as announced. 

A. 5 Calculation for section 03 

For the Franson interferometer of section 0] we must compute R(t,t). Here however, 

S(uj a , LO bl w' a ,Uj' b ) = (l - e^+ a )) (l - e ^+/3)) - e -i«r+a)\ ^ _ g-iKr+flj ^ ^ 
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where we dropped a global factor 4z. Consequently, R2 given in Ij29(l is the sum of 2 6 = 64 integrals, while 
i?4,2 and i?4 ; 3 given respectively by (|35[) and (|36ll are the sum of 2 s = 256 integrals. 

Let's look at i?2- By inspection, one finds that the integrals whose argument is independent of r are four: 
eejllll, that gives a contribution 1; ll|eeee, that also gives a contribution 1; le|eell, whose contribution is 
e l (a+£). and el|llee, whose contribution is e~^ Q+/3 - ) . In these notations, the first two items correspond to the 
products of terms of line two, the last four items correspond to the products within £ . Finally, performing 
the limit AT — > 00 we find 

R 2 (t,t) = 2IJ AB (l + cos(a + f3)) (38) 
that is indeed l|21|l up to a multiplicative factor 2. Immediately then we have also 

i? 41 (r,r) = 4 1 2 J AB (l + cos(a + (39) 
accounting for the first term of the r.h.s. of Ij22(l . 

Moving to Ri t 2, by inspection, one can verify that the only integrals that will not average to zero are those 
associated to the following four products: llll|eeee and eeee|llll, both giving 1; lell|eell, that gives 
e «(a+/3). anc j elee|llee, that gives e~ l ^ a+/3 \ As before, in these notations the first four items represent the 
terms of line two of l|35|) . the last four items correspond to the products within E . After the usual limit, one 
finds 

#4,2 (r, t) +c.c = 2 1 2 J 4 (1 + cos(a + /3)) (40) 
accounting for the second term on the r.h.s. of (|22|l . 

As for ^4,3, again only four integrals out of 256 will not average to zero, namely those associated to lllljeeee, 
eeeejllll, lele|elel and elel|lele; however here, all these contributions give simply 1, so finally 

i?4,3(r,r) - 4I 2 J A J B (41) 

that is the last term in the r.h.s. of (|22Jl . This concludes our demonstration. 
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